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SECTION  1 
INTRODUCTION 


During  the  past  year  or  two,  we  have  attempted  to  make  calcula¬ 
tions  in  which  uncertainties  in  specific  chemical  reaction  rates  could  be 
correlated  to  variations  in  UHF  satellite  communications  degradation.  Ef¬ 
forts  to  this  effect  using  the  ROSCOE  nuclear  weapons  effects  code  were  not 
successful  and  our  attention  turned  to  the  MELT  code.  The  MELT  code  was 
designed  to  describe  the  phenomenology  following  a  high  altitude  nuclear 
burst  out  to  many  hours,  something  the  ROSCOE  code  was  not  designed  to  do. 


The  degradation  of  UHF  and  higher  frequency  signals  at  high  alti¬ 
tudes  is  in  many  ways  qualitatively  different  from  the  types  of  degradation 
encountered  at  lower  frequencies.  Dispersion,  diffraction,  scintillation, 
and  phase  effects  play  a  relatively  more  important  role  at  high  altitudes 
than  at  lower  altitudes  where  collision  induced  attenuation  is  the  dominant 
effect.  At  UHF  and  higher  frequencies,  the  structure  of  the  ionized  region 
causing  the  degradation  becomes  an  important  factor.  Large  scale  regions 
of  ionization  are  known  to  striate  at  high  altitudes  and  the  spacing,  am¬ 
plitude,  and  growth  rate  of  these  striations  are  important  parameters  in 
calculating  the  effect  of  these  regions  of  ionization  on  UHF  and  higher 
frequency  signals  propagating  through  them. 


Theoretically,  a  variation  in  the  chemistry  parameters  could  im¬ 
pact  on  the  growth  and  development  of  striations  as  well  as  on  the  overall 
electron  density,  thus  having  a  compound  effect  on  communications  degrada¬ 
tion.  Unfortunately,  the  striation  parameters  in  present  high  altitude 
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phenomenology  codes  are  specified  in  an  empirical  manner  and  are  not  de¬ 
rived  from  other  physical  parameters.  Under  these  conditions,  a  variation 
in  chemistry  would  not  impact  the  striation  parameters  and  we  had  to  con¬ 
tent  ourselves  with  just  the  effect  of  chemistry  uncertainties  on  either 
the  total  electron  content  within  the  fireball  or  plume  or  the  integrated 
electron  content  along  a  propagation  path. 

The  reactions  and  rate  coefficients  which  are  in  the  MELT  code 
are  treated  in  Section  2.  The  relative  importance  of  these  reactions  is 
discussed  and  the  rate  coefficient  data  are  examined  and  compared  with  the 
coefficients  used  in  MELT.  The  uncertainties  in  these  rate  coefficients 
are  estimated. 

Section  3  deals  with  the  gross  results  of  the  MELT  code  in  the 
area  of  the  relative  importance  of  specific  reactions  and  a  possibly  impor¬ 
tant  deficiency  in  the  code  is  noted. 

Section  4  deals  with  effect  of  the  revised  rate  coefficients  and 
their  uncertainties  on  the  plume  electron  content  as  a  function  of  time. 

Section  5  deals  with  variations  due  to  the  chemistry  uncertain¬ 
ties  of  the  tube  electron  content  (TEC)  along  two  hypothetical  propagation 
paths  and  the  effects  this  might  have  on  the  degradation  of  UHF  and  higher 
frequency  communications.  In  order  to  appraise  possible  scintillation  ef¬ 
fects,  fully  developed  striations  were  assumed. 

In  Section  6  we  give  our  recommendations  for  the  rate  coefficients 
of  the  critical  reactions  involved  in  the  high  altitude  deionization  pro¬ 
cess.  Reactions  which  are  not  now  in  MELT  but  possibly  should  be  included 
are  also  discussed. 
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The  last  section  deals  with  our  recommendation  for  further  work 
in  the  general  area  of  high  altitude  nuclear  induced  chemistry. 


SECTION  2 

REACTIONS  AND  RATE  COEFFICIENTS 


First  we  have  the  electron-ion  recombination  reactions 


The  first  two  reactions  actually  represent  a  complex  series  of  processes 

by  which  an  electron  recombines,  either  by  a  collisional  or  radiative 

transition  into  a  highly  excited  neutral  atom  state  which  is  then  de- 

excited  or  reionized  by  radiative  or  collisional  transitions.  At  low 

electron  densities,  radiative  processes  dominate  and  the  reaction  is 

second  order  in  charged  particle  density.  At  high  electron  densities, 

collisional  processes  dominate  and  the  overall  reaction  is  third  order. 

The  values  used  for  the  rate  coefficient  in  the  MELT  code  for  N+  and 

0+  are  shown  in  Table  1,  where  n  is  the  electron  density  and  T  is 

e  e 

the  electron  temperature  in  Kelvin.  The  rate  coefficients  are  expressed 
as  second  order  in  units  of  cm'Vsec  even  when  collisions  dominate.  These 
values  are  theoretical1  and  very  little  supportive  experimental  data  is 
available.  The  uncertainty,  however,  in  these  values  is  probably  fairly 
small,  particularly  when  the  electron  density  is  low. 
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Table  1.  Radiati ve-Coll i sional  Rate  Coefficients  (in  cm  /sec) 


The  various  dissociative  recombination  reactions,  which  is  shown 
for  NO  in  reaction  3,  have  rate  coefficients  which  are  many  orders  of 
magnitude  faster  than  the  radiative-collisional  rate  coefficients.  The 
rate  coefficie.it  used  for  reaction  3  in  the  MELT  code  is 


Because  of  this  large  difference  in  the  atomic  and  molecular  recombination 
rates,  the  molecular  ions  formed  initially  disappear  almost  immediately 
relative  to  the  atomic  ions  si  ch  as  \'+  and  0+.  The  long-term  disap¬ 
pearance  of  the  atomic  ions  is  controlled  by  either  radiative-col lisonal 
recombination,  when  the  electron  density  is  high,  or  reactions  which  con¬ 
vert  the  atomic  ions  to  molecular  ions  which  then  disappear  almost  immedi¬ 
ately  relative  to  the  remaining  atomic  ions.  Thus,  the  rate  controlling 
step,  when  the  electron  density  is  moderate  to  low,  is  the  conversion  of 
atomic  ions  to  molecular  ions  and  the  values  of  the  dissociative  recom¬ 
bination  rate  coefficients  are  not  critical.  One  of  the  reasons  reaction 
3,  however,  is  important  is  due  to  the  reverse  reaction,  i.e., 


cm' /sec 


n\l 
e  \  e 

500K 

1000K 

2000K 

4000K 

6000K 

107  cm'3 

4.1  E-12 

1.6  E-12 

7.6  E-l 3 

3.8  E-l 3 

— 

2.5  E-l 3 

108  cm'3 

9.3  E-12 

2.6  E-12 

9.7  E-l 3 

4.3  E-l 3 

2.7  E-l 3 

TO9  cm'3 

2.7  E-l 1 

4.8  E-12 

1  .4  E-12 

5.1  E-l 3 

3.1  E-l 3 

1010  cm'3 

1.1  E-l 0 

1.2  E-l 1 

2.3  E-12 

6.7  E-l 3 

3.6  E-l 3 

m11  "3 

1 0  cm 

1.0  E-9 

6.1  E-l  1 

7.1  E-12 

1.5  E-12 

4.0  E-l 3 

m12  -3 

1 0  cm 

9.0  E-9 

3.6  E-l 0 

2.5  E-l 1 

3.1  E-12 

4.0  E-l 3 

(5) 


N  +  0  -*•  N0+  +  e 

The  rate  coefficient  for  this  reaction  carried  in  the  MELT  code  and  ob¬ 
tained  by  detailed  balance  is 

k5  =  3.6  x  10" 12  exp  (-31800/Tn)  cm3/sec  (6) 

where  T^  is  the  neutral  particle  temperature  in  Kelvin.  The  pre¬ 
exponential  factor  is  actually  a  function  of  temperature  and  the  value 
shown  is  for  about  5000  K.  At  the  elevated  temperatures  found  in  the 
fireball  or  plume,  reaction  5  can  be  an  important  source  of  ionization. 

In  addition,  reaction  3  is  known  to  produce  the  excited  state  of  atomic 
2 

nitrogen,  N(“D),  as  a  product.  The  presently  accepted  value  for  the 

2 

fraction  of  the  nitrogen  produced  as  N(  D)  is  0.75.  If  we  consider 
2 

N(  D)  only  in  reaction  5,  i.e., 

N(2D)  +  0  -+  N0+  +  e  (7) 

we  obtain  a  rate  constant  of 

ky  =  3.0  x  io"13  C^qO  1/2  exp  (-4300/T)  cm3/sec  (8) 

This  can  produce  an  even  greater  amount  of  ionization  at  even  lower  tem- 

4 

peratures  than  reaction  5  using  N  in  the  S  ground  state.  To  get  an 

2 

appreciable  effect,  however,  the  concentration  of  N(  D)  must  be  signifi¬ 
cantly  larger  than  that  at  thermal  equilibrium.  If  N(“D)  is  in  thermal 
equilibrium  with  the  ground  state,  then  its  effect  is  not  significant  and 
is  included  in  the  rate  for  reaction  5  when  the  total  nitrogen  atom  con¬ 
centration  is  used. 
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N(  D)  is  produced  in  nonequilibrium  amounts  in  and  near  the  fire¬ 
ball  primarily  from  dissociative  recombination  of  N*  and  N0+,  i.e., 


N*  +  e  -+  N(2D)  +  N(4S) 


N0+  +  e  ->  C .  75  N(2D)  +  0.25  N(4S)  +  0 


(9) 

(10) 


2  4 

The  proportions  of  N(  D)  produced  relative  to  the  ground  state,  N(  S) , 

2 

are  somewhat  uncertain  but  probably  not  by  more  than  20  percent.  N(  D)  is 
probably  produced  also  by  the  reaction 


n+  +  o2  -*  n(2d)  +  o; 


(ID 


though  no  experimental  evidence  is  available.  This  is  likely  because  the 
above  reaction  is  very  nearly  resonant  and  charge  exchange  reactions  prefer 
channels  which  are  resonant.  Even  in  regions  where  the  0^  is  ultimately 
destroyed,  this  reaction  should  contribute  at  early  times. 

2 

Whether  these  large  nonequilibrium  concentrations  of  N(  D)  re¬ 
main  over  long  periods  of  time  depends  in  large  measure  on  whether  molecular 

2 

oxygen  is  present  in  sufficient  quantities  to  quench  the  N(  D)  by  the 
reaction 


The  NO 


N(2D)  + 
produced  will 
N(2D)  + 


02  -*  NO  +  0  . 

help  deplete  the 
NO  -+  N2  +  0  . 


N(2D) 


by  the  reaction 


- 12 

The  rate  constants  for  these  reactions  are  7  *  10  and  6  x  10 
respectively. 


(12) 


(13) 

3. 

cm  /sec. 
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In  regions  where  the  0^  has  been  destroyed  or  is  not  present 

2 

in  sufficient  quantity  the  quenching  of  N(  D)  proceeds  more  slowly.  The 

quenching  by  electrons  is  fast  when  the  electron  density  is  large  but  the 

electron  concentration  may  not  remain  large  enough  for  enough  time  to 

2 

significantly  affect  the  N(  D)  depletion.  The  quenching  rate  for  molecular 
nitrogen  is  small.2  The  rate  for  atomic  oxygen  is  somewhat  uncertain  but 
is  probably  of  the  order  of  magnitude  of  10  “  cm  /sec.3  Thus,  in  the 
region  of  150  to  250  kilometers  altitude  the  characteristic  decay  time  for 
N(2D)  will  be  of  the  order  of  100  to  1000  seconds.  If  the  N(“D)  is 
several  orders  of  magnitude  above  its  equilibrium  value  it  may  take  an 
order  of  magnitude  longer  than  its  characteristic  decay  time  to  reach 
equilibrium.  Thus,  under  the  right  conditions,  reaction  7  can  be  a  signi¬ 
ficant  source  of  ionization.  This  was  demonstrated  when  the  results  of 
the  high  altitude  chemistry  model  written  for  ROSCOE  were  compared  to  a 
number  of  benchmark  calculations.  In  one  case,  that  of  a  parcel  of  air 
at  about  150  km  altitude  close  to,  but  not  in,  the  fireball,  the  ionization 
obtained  over  a  long  period  of  time  by  the  ROSCOE  model  was  significantly 

greater  than  that  any  of  the  other  codes  had  obtained.  This  was  due  to 
2 

the  N(  D) ,  which  was  carried  in  the  ROSCOE  model  and  allowed  to  deviate 
from  equilibrium. 

The  experience  with  the  ROSCOE  model,  however,  indicated  that 

this  effect  was  not  widespread  but  limited,  at  least  during  the  first  10 

2 

minutes,  to  regions  near  the  burst  in  which  large  amounts  of  N(  D)  are 
produced  and  the  ambient  0£  is  destroyed  by  the  energy  input  from  the 
weapon.  At  altitudes  above  200  or  250  kilometers  and  in  the  fireball, 
the  0+  dominated  the  ionization  to  the  extent  that  the  N0+  was  not 
significant.  A  systematic  study  of  this,  however,  was  not  made. 

2 

Unfortunately,  the  MELT  code  does  not  carry  N(  D)  as  a  separate 
species  and  therefore  we  could  not  assess  the  importance  of  reaction  7  at 
late  times.  In  addition,  MELT  does  not  do  neutral  chemistry  and  at 


12 


80  minutes  all  neutral  species  concentrations  are  replaced  by  ambient 
concentrations . 

Again  one  of  the  main  reasons  for  the  importance  of  reaction  3 
is  the  importance  of  the  reverse  reactions,  the  rates  of  which  are  usually 
obtained  from  the  forward  rates  by  detailed  balance.  Another  reason  for 

~h  +  + 

the  importance  of  NO  recombination,  as  contrasted  with  and  0? 

recombination,  is  due  to  the  fact  that  nuclear  bursts  produce  nitric  oxide. 
This  nitric  oxide,  or  NO,  will  persist  in  the  atmosphere  for  a  long  period 
of  time  after  the  bursts  have  occurred  and  after  the  other  atmospheric 
species  concentrations  have  relaxed  to  ambient  or  near  ambient  values. 

This  NO  can  be  ionized  by  the  Lyman  a  radiation  from  the  sun  and  will 
increase  the  natural  ionization  in  the  atmosphere  for  long  periods  of  time 
after  a  burst.  The  magnitude  of  this  effect  depends  on  the  amount  of  the 
extra  NO  produced,  the  production  rate  of  N0+  by  Lyman  a,  and  the  rate 
of  disappearance  of  the  NO+.  Thus  the  recombination  rate  of  N0+  is  also 
important  for  this  reason. 

The  value  of  the  rate  coefficient  for  N0+  dissociative  recomb¬ 
ination  at  300K  is  fairly  well  established.  At  higher  temperatures  the 
situation  is  confused.  Measurements  in  which  the  overall  temperature  was 
varied  yield  a  T  *  dependence.  Measurements  in  which  only  the  electron 
temperature  was  varied  yield  a  T~‘^  dependence.4  At  1000K  and  3000K 
this  amounts  to  factors  of  2  and  4  respectively  in  the  rate  coefficient. 

The  electron  density  will  depend  at  most  on  the  square  root  of  the  re¬ 
combination  coefficient  and  thus  this  does  not  represent  a  very  serious 
uncertainty,  particularly  if  the  real  values  are  indeed  bracketed  by  the 
two  temperature  dependences.  In  addition  the  measurement  discrepancy 
should  be  resolved  within  the  chemical  kinetics  community  with  a  reason¬ 
able  amount  of  time. 
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We  come  now  to  the  important  exchange  reactions  which  convert 

the  atomic  ions  N+  and  0+  to  molecular  ions  such  as  0*  and  N0+. 

It  has  been  consistently  demonstrated  in  the  past  that  in  most  regions 

these  reactions  are  the  controlling  factor  in  the  deionization  process. 

In  regions  where  0^  is  present  in  sufficient  quantities,  the  N+  will 
+  + 

be  converted  to  NO  and  0^  by  the  reactions 


N+  +  02  -+  0*  +  N(2D) 

(14) 

N+  +  09  -v  N0+  +  0 

(15) 

The  best  value  at  this  time  for  the  sum  of  the  rate  coefficients  of  these 

two  reactions  is  6.1  x  10  cm  /sec.5’6  (Because  the  MELT  code  does  not 

explicitly  carry  the  species  0+,  only  reaction  15  is  carried  with  a  rate 

-10  3  2 

coefficient  of  5  x  10  cm  (sec).)  The  uncertainty  for  the  total  rate 

coefficient  at  300K  is  small  enough  to  be  unimportant  in  the  context  of 

high  altitude  nuclear  weapons  effects  and  the  dependence  on  temperature 

in  the  region  of  interest  is  slight.  The  other  uncertainties  regarding 

these  reactions  are  in  the  branching  ratio  between  reactions  14  and  15 

2 

and  the  amount  of  N(  D)  produced.  An  even  split  in  the  branching  ratio 

is  usually  assumed  though  the  most  recent  measurement5  yielded  51  percent 

and  43  percent  for  reactions  14  and  15,  respectively,  with  the  remaining 

6  percent  going  to  0+  and  NO  as  the  products.  The  main  impact  of  these 

2 

uncertainties  would  be  in  the  N(  D)  and  NO  production.  Since  0,  is 
present,  however,  subsequent  reactions  would  yield  close  to  the  same  result 
regardless  of  the  specific  initial  channel  assumed,  i.e.,  the  production 
of  NO  and  three  oxygen  atoms.  The  uncertainties  are,  therefore,  unimpor¬ 
tant  for  our  purpose. 

In  regions  where  molecular  oxygen  is  not  present  in  sufficient 
quantities  to  accommodate  the  N+,  such  as  at  very  high  altitudes  and  in 
and  near  the  fireball,  the  charge  exchange  reaction 


(16) 


N+  +  0  -*■  0+  +  N 

is  a  necessary  first  step  in  neutralizing  the  charge.  The  0+  formed  can 

then  enter  into  reactions,  to  be  discussed  later,  which  yield  molecular 

ions.  Thus,  if  large  quantities  of  N+  are  formed,  reaction  16  is  of 

critical  importance  in  the  deionization.  Of  the  significant  reactions, 

reaction  16  is  also  the  one  with  the  greatest  uncertainty.  The  MELT  code 

-l9  3 

carries  a  constant  rate  coefficient  rate  of  10  cm  /sec  up  to  a  value 

of  0.2  eV  or  about  2300K  for  the  sum  of  the  ion  and  neutral  temperatures. 

Above  that  the  rate  coefficient  is  taken  to  be  a  linear  function  of  the 

sum  of  the  ion  and  neutral  temperatures.  The  constant  rate  coefficient 

at  low  temperatures  is  a  guess.  The  values  at  high  temperatures  are  based 

on  an  extrapolation  of  cross  section  measurements  at  energies  of  5  eV  and 

above.  More  recent  experiments  have  extended  the  data  down  to  0.5  eV, 

which  is  much  closer  to  the  temperature  regime  of  interest.7  This  data  is 

presented  in  Figure  1  (circles)  and  is  given  in  terms  of  the  measured  cross 

section,  a,  times  the  relative  velocity,  v,  of  the  colliding  particles  in 

the  center  of  mass  system  versus  the  center  of  mass  energy.  The  solid  line 

represents  a  fit  to  these  points  between  0.5  eV  and  3.5  eV  energy.  Though 

3 

the  quantity  cv  has  the  dimensions  of  cm  /sec,  it  is  not,  strictly  speak¬ 
ing,  a  rate  coefficient.  However,  when  the  fit  to  these  values  is  inte¬ 
grated  over  a  Maxwellian  velocity  distribution,  we  obtain  a  rate  coefficient 
equal  to  the  solid  line  for  the  temperature  scale  shown  in  Figure  1. 

This  yields 


k16  =  3  x  10'13  (-!-)  cm3/sec  (17) 

The  triangular  point  in  Figure  1  represents  a  preliminary  upper  limit  for 
the  cross  section  at  0.19  eV.  This  value  was  not  used  in  obtaining  Equa¬ 
tion  17  but  is  shown  in  order  to  indicate  that,  if  true,  it  does  not  re¬ 
sult  in  any  inconsistency. 
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Equation  17  yields  values  for  which  are  within  20  percent 

of  the  values  used  in  MELT  above  800K.  Because  of  the  scatter  in  the  data, 
the  total  uncertainty  above  about  4000K  is  easily  a  factor  of  three,  (vhere 
the  fit  depends  strongly  on  the  extrapolation  of  these  data  values,  the  un¬ 
certainty  increases.  At  1000K,  an  order  of  magnitude  or  more  uncertainty 
must  be  considered.* 

I 

The  0+  ion  is  the  most  important  ion  throughout  most  of  the 
space  above  about  120  km,  particularly  if  we  assume  the  N  can  charge 
exchange  with  atomic  oxygen  quickly  enough.  In  the  absence  of  any  large 
source  of  ionization,  0+  will  be  the  last  ion  to  disappear  and  will  thus 
dominate  the  late  time  ionization.  The  two  primary  reactions  involving 
0+  are 


o+  +  o?  -*  0*  +  0 


0+  +  N?  -*  N0+  +  N 


(18) 
(19) 

The  rate  coefficient  carried  in  the  MELT  code  for  reaction  18  is  2  *  10 
cm'Vsec  for  all  temperatures.  This  is  about  the  measured  value  at 
300K5’ 8’  9  ’ 1 1  ’ 1 4  but  higher  temperature  measurements8’3’10’11’14  show  a 
distinct  decrease  in  the  rate  coefficient  up  to  about  1800K  at  which  point 
the  rate  coefficient  is  less  than  one-half  that  at  300K.  The  rate  coef¬ 
ficient  rises  rapidly  again  above  1800K  and  at  about  3000K  it  is  2  *  10 
again.  At  higher  temperatures,  however,  there  is  not  likely  to  be  any 
07  at  these  altitudes  so  we  need  not  be  strongly  concerned  with  the  data 


11 


Very  recently  an  analysis  of  the  Atmospheric  Explorer  data  has  yielded 
a  value  for  k^  of  2,2  x  10"^  cm-V sec  at  840K  with  a  claimed  accu¬ 
racy  of  ±  20  percent.17  This  value  has  been  added  to  Figure  1.  The 
dashed  line  represents  a  low  temperature  fit  discussed  in  a  footnote  in 
Section  6. 
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above  3000K.  The  uncertainty  is  fairly  small  at  300K  (20  percent)  but 
because  of  scatter  in  the  data,  increases  to  about  a  factor  of  two  at 
about  1000K. 8,9  0 jl  4  This  is  a  tolerable  uncertainty  if  the  alternate 

path  for  0  depletion,  reaction  19,  is  strongly  dominant.  If,  however, 
reaction  18  makes  up  a  significant,  even  if  not  dominant,  fraction  of  the 
0  depletion,  this  uncertainty  becomes  significant. 

Reaction  19  is  the  most  important  reaction  involved  in  high  alti¬ 
tude  deionization.  The  rate  coefficient  in  the  MELT  code  is,  in  units  of 
cm  /sec. 


19-1.2*10-12 

T*  <  300K 

19  ’  >-2  *  10"12 

300K  <  T*  <  1200K 

19  ■  '-2  *  10-12( 

T*  \2 

1200K<  T*<  10000K 

1200  / 

ig  =  8.26  x  10'11 

i  T*  \ 

T*  >  10000K 

\ioooo/ 

where  T*  is  the  greater  of  the  ion  temperature  and  the  N2  vibrational 
temperature.  This  points  up  not  only  the  strong  dependence  on  translation 
temperature  above  1200K  but  also  the  strong  dependence  on  the  vibrational 
state  of  the  reacting  nitrogen  molecule.  Figures  2l?  and  38  show  the 
experimental  data  for  as  a  function  of  N2  vibrational  temperature 

and  translational  energy.  This  is  probably  the  best  data  available.  The 
constant  value  of  the  rate  coefficient  up  to  about  a  thousand  degrees 
nitrogen  vibrational  temperature  (Figure  2)  is  attributed  to  the  fact 
that  the  rate  coefficient,  at  300K  translational  temperature,  of  0+  with 
nitrogen  in  the  v  =  0  and  v  =  1  vibrational  states  is  the  same  (1.2 
X  10  cm  /sec)  and  the  v  =  2  state  has  not  been  appreciably  populated 
yet.  The  rapid  increase  in  the  rate  constant  above  1000K,  when  the  v  >  1 
levels  progressively  become  populated,  is  taken  to  mean  that  the  rate 


18 


)  50  ’00  40  9C  200  2*0  230 

£,'N  tTd)  A<gon  Buffer 


1 

i 


,a 

8 

a 

u 

3  m 

€ 
|S?  ? 
?> 
J 

I®  pi 

7 
>  U>| 


1TTT  T  T  V  } - T" 


l 


if  I  S!g  S  8sSe  I 

CD  ®  -  oj  ®<\i*o*en  ? 

®  cC5c3E3cS<3rt| 

C\J  o>  •  •  J  ■  *  •  •  ° 


Z 

♦ 

b 


z 

» 

o 


J 


.1 

1 1 1 1  i  i  i  i _ !i  l_i  j i i i i _ JjLU-U-LJa 

“  V  T  ?  T  Trt 

O  O  O  O  O 


|3  " 

E 

u 
UJ 


fa  fa 

X  X 


—  */■»**■* 


(tS  1 91039101 U  ?UJ0)  IUO1SUO0  9|0H 


C  O 
O  *r“ 
4-> 

4- >  CD 
U  C 

fT3  *r- 
<U 

5- 

c 

4-  O 
O  «r- 


C  O 
0) 

r-  C 

cj  o 


CD 
i- 
I  3 
ra  +-> 

S-  03 
CD  S- 
Q_  CD 

E  cl 

cu  e 

4->  <L> 
4-> 

<TJ  CJ 
C  *r- 

O  +-> 
•r~  O) 
-M  C 
H3  *r- 

i.  ^  « 
_Q 


4-  O 
CD  C 
O  3 
<J  4- 

<D  n3 

rO  to 
C£  fT3 


CD 

3 

CD 


C\f*£ 
Z  CO 
—  03 
C\J 

>3  II 
CD 

CD 
CD  S- 
C  3 
CD  -«-> 


O 

CO 

r— 


f— 

CD 

(V 

S- 

C3  C 

3 

r—  O 

4-> 

•r— 

03 

C  4-> 

S- 

O  fC3 

CD 

•r-  S- 

Ql 

4->  -Q 

E 

(J  *i- 

CD 

03  > 
CD 

4-> 

S- 

o 

CNi 

•r— 

4-  Z 

-»-> 

O 

CD 

C 

4J  4- 

•r- 

C  o 

CD 

£ 


o  o 

•r-  *r- 

4-  -♦-> 
4-  O 
CD  C 
O  p 
U  4- 

CD  <T3 
-♦-> 

ro  to 
Of  ro 


CD 

U 

3 

CD 


(D 

S- 

3 

4-> 

fO  * 
U  ^ 
CD  O 
CLO 

E  ro 
CD 

♦->  II 


19 


constant  with  N~  in  the  v  =  2  and  higher  levels  is  much  larger  than 

that  for  v  =  0  and  1.  Analysis  of  the  data  in  Figure  2  yields  a  value 

5  x  10  11  cm'Vsec  for  the  reaction  with  the  N?  in  the  v  =  2  state.12 

The  rate  constants  continue  to  increase  with  increasing  values  of  v  and 

level  off  at  about  v  =  6  at  a  value  of  about  5  x  10  0  cm'Vsec.12  It 

should  be  emphasized  that  the  rate  coefficients  in  Figure  2  apply  only  when 

the  translational  temperature  is  300K.  When  the  ^  is  in  the  lowest 

vibration  state  and  the  kinetic  temperature  is  varied  we  have  the  situation 

in  Figure  3.  We  see  that  the  value  of  the  rate  constant  decreases  from 

its  value  at  300K  with  increasing  temperature  (T  =  7730KEcm)  and  reaches 

a  minimum  value  of  about  5  x  10  ^  cm'Vsec  at  about  KE  =  0.16  eV  or 

cm 

1200K.  In  this  temperature  range,  when  the  higher  vibrational  levels  (v  >  1) 
of  nitrogen  are  not  significantly  excited,  it  is  the  translational  temper¬ 
ature  which  controls  the  variation  of  the  rate  coefficient.  Thus,  the  use 
of  a  constant  rate  coefficient  between  300K  and  1200K  is  wrong  by  up  to  a 
factor  of  2.4.  This  can  be  significant  since  the  ambient  temperature  of 
the  atmosphere  encompasses  just  this  range  above  100  km  altitude.  Above 
200  km  altitude,  the  range  of  ambient  temperature  is  600K  to  1200K.  Thus, 
the  rate  coefficient  of  reaction  19  used  in  obtaining  most  of  the  late 
time  MELT  results  is  wrong  by  about  a  factor  of  two.  We  will  discuss  this 
more  later. 

Above  about  2000K  the  rate  coefficients  as  a  function  of  nitrogen 
vibrational  temperature  and  kinetic  temperature  are  approximately  the  same. 
However,  this  does  not  mean  these  values  are  correct  to  use  when  both  the 
ion  and  vibrational  temperatures  are  in  this  range.  To  repeat,  the  data 
in  Figure  2  is  for  a  kinetic  temperature  of  300K  only  and  the  data  in  Fig¬ 
ure  3  is  only  for  the  lowest  vibration  state.  The  rate  coefficients  when 
both  temperatures  are  high  may  be  quite  different.  There  is,  unfortunately, 
no  experimental  data  in  this  regime  to  guide  us. 
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Theoretical  calculations  have  been  made  for  the  dependence  of 
the  rate  coefficient  as  a  function  of  both  kinetic  and  vibrational  temper¬ 
atures  up  to  7000K . 1 3  The  results  agree  fairly  well  with  the  data  in 
Figure  2,  when  only  the  vibrational  temperature  is  varied,  but  not  too 
well  with  the  data  in  Figure  3,  when  only  the  translational  temperature 
is  varied.  The  calculations  are  probably  a  good  guide  to  the  qualitative 
behavior  of  the  rate  coefficient  when  both  temperatures  are  varied  but  a 
fairly  significant  uncertainty  still  exists  quantitatively. 

Another  source  of  uncertainty  is  apparent  in  Figure  3  and  applies 
to  other  rate  coefficient  measurements  as  well.  The  measured  coefficients 
are  different  for  different  buffer  gases  (helium  and  argon)  used  in  this 
particular  set  of  experiments.  Extrapolation  of  the  argon  data  to  lower 
energies  shows  significant  differences  down  to  about  600K.8  The  helium 
data  more  closely  approximates  the  case  for  a  Maxwellian  speed  distribu¬ 
tion.  However,  the  procedure  of  relating  the  ion  speed  distribution  brought 
about  by  the  application  of  an  electric  field  to  an  effective  thermal  tem¬ 
perature  introduces  uncertainties,  particularly  when  the  internal  excitation 
of  the  reactants  plays  a  role  in  the  reaction.8 

Another  problem  is  how  all  of  the  rate  coefficient  data  is  applied 
to  a  case  in  which  the  translational  temperatures  of  the  reactants  are  dif¬ 
ferent.  This  is  a  different  problem  than  when  the  translational  and  vi¬ 
brational  temperatures  are  different.  Analysis15  indicates  that  one  may 
use  a  single  translational  temperature  which  is  the  mean  of  the  transla¬ 
tional  temperatures  of  the  reactants,  each  of  which  is  weighted  by  the 
inverse  of  the  molecular  weight.  For  the  0+  +  N_  reaction,  this  yields 
a  T  which  is  equal  to  0.64  T.  +  0.36  T  .  For  the  0+  +  0,  reaction, 
we  have  T  =  0.67  T.  +  0.33  Tn>  for  the  N+  +  C>2  reaction,  T  =  0.70  T. 

+  0.30  T  and  for  the  N+  +  0  reaction,  T  =  0.53  T.  +  0.47  T  . 

n  in 
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SECTION  3 

MELT  CHEMISTRY  RESULTS 


We  shall  discuss  in  this  section  two  separate  calculations  using 
the  MELT  code.  In  the  first  case  the  MELT  code  was  run  for  a  megaton  range 
burst  at  200  km  altitude  to  a  time  of  about  30  minutes.  The  procedure 
beyond  that  time  was  to  choose  the  central  striation  in  the  plume  and  carry 
it  forward  to  about  12  hours  using  the  MELT  chemistry  routines  in  conjunc¬ 
tion  with  crude  hydrodynamics.18  Figure  4  contains  the  results  of  this  cal¬ 
culation  presented  in  a  way  which  shows  the  relative  importance  of  the 
different  reactions  we  have  discussed.  This  figure  shows  the  number  of 
reaction  events,  on  a  relative  basis,  which  have  taken  place  up  to  a 
specific  time.  Since  the  values  given  at  any  time  are  proportional  to 
the  integrated  number  of  reaction  events  for  each  reaction,  the  slope  of 
the  curves  at  any  given  time  is  indicative  of  the  relative  importance  of 
each  reaction  at  that  particular  time.  We  see  that  this  figure  tends  to 
contradict  the  assertions  we  made  in  the  previous  section  concerning  which 
reactions  are  the  most  important.  The  main  bulk  of  the  ionization  is 
removed  by  the  radiative  recombination  reactions  (a  reaction  which  results 
in  the  formation  of  a  molecular  ion  is  considered  to  be  an  ion  removal 
reaction  because  of  the  rapid  subsequent  dissociative  recombination)  and 
up  to  several  hours  these  reactions  are  the  fastest  ion  removal  reactions. 
The  0+  +  N0  reaction  (reaction  19)  is  never  dominant  and  contributes  no 
more  than  the  0+  +  0^  reaction  (reaction  18).  This,  however,  is  not  a 
typical  case.  This  calculation  dealt  only  with  the  central  tube  or  stri¬ 
ation  of  the  plume  for  a  burst  in  which  the  total  electron  content  was 
initially  much  higher  than  is  usual.  Since  the  radiative  recombination 
rate  depends  on  the  square  of  the  ion  or  electron  density  and  the  other 


Figure  4.  Relative  number  of  reactions  (compared  to  the  total  number  of 
deionizing  reactions  during  the  12  hour  simulation)  of  the  va 
chemical  reactions. 


reactions  depend  only  on  the  first  power  of  the  ion  density,  the  radiative 
recombination  is  enhanced  relative  to  the  ion-neutral  reactions  when  the 
ionization  level  is  high.  This,  however,  is  only  part  of  the  reason. 
Another  reason  is  that  because  of  a  complicated  interplay  of  the  hydro¬ 
dynamics,  plasma-neutral  interaction  and  chemistry  during  the  time  span 

16  _  2 

when  the  tube  electron  content  is  high  (>  2  x  10  cm  ),  the  bulk  of 
the  electron  content  is  held  at  high  altitudes.  At  these  altitudes  there 
are  ve  few  molecules  so  the  radiative  recombination  tends  to  be  dominant 
in  these  regions. 

The  second  calculation19  gives  results  which  we  believe  are  more 
typical  of  the  overall  plume  chemistry.  This  calculation  was  a  full  blown 
MELT  calculation  for  a  lower  yield  burst  in  the  megaton  range  at  250  km 
altitude  which  was  run  to  ten  hours.  Figure  5  shows  information  from  this 
run  similar  to  that  in  Figure  4  except  that  the  total  number  of  reactions 
are  plotted  and  the  integration  is  over  the  whole  plume  rather  than  over 
just  the  central  striation.  We  see  that  the  primary  ion  removal  reaction 
at  most  times  is  reaction  19,  which  accounts  for  about  46  percent  of  all 
the  ions  removed  at  ten  hours.  In  addition,  the  bulk  of  the  N+  ions  are 
depleted  by  reaction  16.  This  is  in  agreement  with  what  past  experience 
has  indicated  would  be  the  important  reactions. 

Figure  5  shows  certain  discontinuities  at  about  80  minutes,  par¬ 
ticularly  in  the  reactions  involving  0r  This  is  due  to  a  crude  artifice 
in  this  particular  MELT  calculation ,  the  purpose  of  which  is  to  compensate 
for  the  lack  of  neutral  diffusion  in  the  code.  At  about  80  minutes,  it  is 
assumed  that  diffusion,  were  it  in  the  code,  would  have  brought  the  atmo¬ 
sphere  back  to  near  normal  with  respect  to  neutral  species  concentrations 
and  so  all  neutral  concentrations  are  replaced  by  their  corresponding  am¬ 
bient  concentrations.  Thus,  at  80  minutes,  the  regions  in  the  plume  which 
have  been  depleted  of  ambient  molecular  oxygen,  for  instance,  suddenly 
have  it  replaced.  There  is  at  least  one  possible  objection  to  this 
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Number  of  Reactions  (Units  of  10 


Figure  5.  Chemical  reactions  as  a  function  of  time  after  burst 
(2.5  minutes  to  10  hours). 
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procedure.  We  have  discovered  from  calculations  regarding  low  altitude 
fireballs  (<100  km),  in  which  entrainment  of  surrounding  ambient  air  was 
included,  that  the  initial  molecular  oxygen  entrained  was  dissociated  by 
the  atomic  nitrogen  in  the  fireball.  The  reactions  accomplishing  this  were 

N  +  09  -*•  NO  +  0  (21) 

N  +  NO  -*  N9  +  0  .  (22) 

It  is  not  until  the  atomic  nitrogen  is  used  up  in  this  or  other  ways  that 
the  molecular  oxygen  concentration  can  build  up  to  appreciable  values. 

The  situation  in  the  MELT  case  is  somewhat  different  as  to  temperature 
and  density  but  the  plume  does  contain  atomic  nitrogen  and  a  certain 
amount  of  0-,  diffusing  into  the  plume  would  be  destroyed  and  the  process 
of  the  0.,  replacement  would  be  prolonged.  This  could  eliminate,  or  at 
least  mitigate,  the  enhancement  of  the  reactions  involving  0 2  seen  in 
Figure  4.  Since  this  effect  would  increase  the  amount  of  atomic  oxygen 
we  might  expect  instead  an  enhancement  of  reaction  16  (N+  +  0  -*•  0+  +  N) . 

This  process  would  also  increase  the  NO  concentrations  in  the  plume, 
the  possible  effects  of  which  were  mentioned  previously.  Since  the  MELT 
code  does  not  carry  any  neutral  reactions  as  well  as  not  carrying  diffusion, 
it  is  difficult  to  determine  how  important  this  effect  is. 

Because  of  the  assumption  of  ambient  O2  concentrations  after 
80  minutes,  the  reactions  of  0+  and  N+  with  0o  (reactions  18,  14, 
and  15)  become  significant.  Reaction  IS,  in  fact,  becomes  dominant  after 
about  eight  hours.  It  is  difficult  to  say  how  significant  this  result  is 
considering  the  discussion  of  the  last  paragraph  and  the  fact  that  these 
reactions  also  tend  to  deplete  the  02< 
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SECTION  4 


EFFECT  OF  NEW  RATE  COEFFICIENTS  AND  UNCERTAINTIES 


In  Figure  6  are  plotted  the  MELT  results  for  the  plume  total  elec¬ 
tron  content  (solid  line)  as  a  function  of  time.  The  altitudes  and  temper¬ 
atures  at  which  most  of  the  reactions  depleting  the  electrons  are  taking 
place  are  as  yet  largely  unknown  although  further  analysis  and  examination 
of  the  MELT  outputs  will  eventually  yield  this  information.  The  exponential 
decay  time  of  1.7  hours,  or  6000  seconds,  would  indicate  an  altitude  region 
of  under  400  km.  This  is  not  surprising  since  this  is  where  the  bulk  of  the 
neutral  species  is  located.  The  temperatures  (the  electron,  ion,  neutral 
and  N2  vibrational  temperatures  are  all  likely  to  be  different)  at  which 
the  reactions  are  taking  place  are  not  easily  estimated  without  further 
analysis  of  the  MELT  outputs.  The  fall  of  the  plasma,  after  the  initial 

rise,  into  the  neutrals  causes  a  renewed  heating  which  may  yield  kinetic  ion 

4 

temperatures  as  high  at  10  K.  Actually,  the  plasma  rises  and  falls  twice 
before  stabilizing  at  about  1  1/2  or  2  hours.  After  this  time  the  temper¬ 
atures  are  probably  fairly  close  to  those  assumed  for  the  ambient  ionosphere. 
These  were  assumed  to  be  1000K  for  the  heavy  particle  kinetic  temperature 
and  1600K  for  the  nitrogen  vibrational  temperature  for  altitues  above  about 
200  km. 

To  get  some  idea  of  what  might  be  the  effect  of  the  different 
rate  coefficients  and  uncertainties  we  have  discussed  we  have  used  what 
we  considered  the  best  values  for  the  reaction  rate  coefficients,  and 
plotted  in  Figure  6  the  electron  content  of  the  plume  versus  time  using 
these  values.  We  have  assumed  that  the  changes  in  the  chemistry  do  not 
affect  the  plasma  and  hydrodynamic  motions  and  we  believe  this  is  a  good 
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assumption  as  long  as  we  restrict  ourselves  to  times  after  the  plume  has 
stabilized.  We  have  taken  this  stabilization  time  to  be  100  minutes  and 
have  applied  the  new  reaction  rates  after  this  time.  This  was  a  simple 
procedure  since  the  total  electron  exponential  decay  time  was  fairly  constant. 
As  long  as  the  ion-molecule  reactions  remain  dominant  we  would  expect  the 
decay  time  to  remain  constant  in  time  even  though  its  value  changes  with 
changes  in  the  rate  coefficients. 

With  the  1600K  nitrogen  vibrational  temperature  the  rate  coef¬ 
ficient  in  MELT  for  reaction  19,  the  reaction  of  0  with  N.,,  is  2.1  * 

10  cm3/sec.  The  reaction  we  have  used  is  based  on  the  fact  that  at 
1600K  about  1.4  percent  of  the  N7  is  in  the  v  =  2  vibrational  state, 

0.2  percent  in  the  v  =  3  state,  and  the  remainder  in  the  V  =  1  and 
v  =  2  states.  From  Figure  2,  we  have  the  rate  constant  for  N-  in  the 
v  =  1  and  v  =  2  states,  at  1000K  translational  temperature,  is 
5  x  10  *3  cm3/sec.  We  assume  that  the  rate  constant  for  v  =  2  at  1000K 
is  the  same  as  at  300K,  i.e.,  is  equal  to  5  *  10  ^  cm3/sec.  Since  1.4 

percent  of  the  N_  is  in  the  v  =  2  state,  we  have  a  contribution  to  the 
1  - 13  3 

overall  coefficient  of  7  x  10  cm  /sec.  Similarly,  the  v  =  3  state 
contributes  2  x  10* 12  cm3/sec.  The  total  value  is  1.4  k  10  *'  cm3/ sec 
or  two-thirds  of  the  rate  coefficient  used  in  MELT. 


The  value  we  have  used  for  the  rate  coefficient  of  reaction  18, 
i.e.,  the  0  +  02  reaction,  at  1000K  is  10"11  cm3/sec.  This  is  half 
the  value  used  presently  in  the  MELT  code. 

At  1000K  the  rate  coefficient  we  use  for  reaction  16,  i.e.,  the 
N  +0  charge  exchange  reaction,  is  the  same  as  the  value  used  in  MELT. 
Although  this  reaction  does  not  directly  produce  an  ion-electron  recombina¬ 
tion  nor  a  molecular  ion  which  can  immediately  recombine,  it  does  produce 
much  of  the  0  which  then  feeds  into  reactions  18  and  19.  A  severe  decrease 
in  the  rate  coefficient  of  reaction  16  could  cause  this  reaction  to  control, 
at  some  point,  the  overall  rate  of  ion  decay.  An  examination  of  Figure  5. 
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shows,  however,  that  the  supply  of  N  is  essentially  depleted  by  about 
3  to  4  hours.  A  severe  decrease  of  a  factor  of  two  or  more  could  have  a 
significant  effect  and  this  is  a  possibility  considering  the  uncertainty 
in  this  rate  coefficient.  A  substantial  increase  in  the  rate  coefficient 
of  reaction  16,  on  the  other  hand,  would  have  little  overall  effect  in  this 
case  since  reactions  18  and  19  are  controlling  the  0+  ion  decay. 

The  rate  coefficients  for  the  other  reactions  involved  were  not 
changed  from  the  values  in  MELT.  The  changes  in  reactions  18  and  19  were 
enough,  however,  to  change  the  exponential  decay  time  from  1.7  hours  to 
2.2  hours.  The  effect  of  this  on  the  total  plume  electron  density  is  shown 
by  the  dashed  line  in  Figure  6.  At  7  hours  we  have  a  difference  of  about 
a  factor  of  two  and  at  10  hours  a  difference  of  almost  a  factor  of  three. 

The  effect  of  the  uncertainties  in  these  rate  coefficients  can 
be  estimated  in  a  similar  manner.  We  have  already  pointed  out  that  the 
differences  and  scatter  in  the  data  from  different  experiments  yields  about 
a  factor  of  two  uncertainty  in  the  rate  coefficients  as  a  function  of  kinetic 
temperature  for  reactions  18  and  19.  The  lack  of  experimental  data  for  re¬ 
action  19  concerning  the  variation  of  the  rate  coefficient  with  kinetic 
temperature  for  the  higher  vibrational  states  certainly  causes  at  least 

an  uncertainty  of  a  factor  of  two  in  this  contribution  to  the  coefficient. 

To  these  uncertainties  must  be  added  the  uncertainty  in  reaction  16, 
the  N+  +  0  charge  exchange  reaction,  which  affects  the  ion  decay  time  in 
an  indirect  way.  As  we  mentioned  before,  the  effect  of  this  uncertainty 
will  be  more  pronounced  when  the  coefficient  is  decreased  than  when  it  is 
increased.  A  factor  of  two  uncertainty  at  about  1000K  either  way  is  not 
unreasonable  for  this  difficult  to  measure  rate  coefficient. 

The  effect  of  these  uncertainties  on  the  plume  electron  content 
is  shown  by  the  dotted  lines  in  Figure  6.  The  electron  content  due  to 
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increasing  the  rates  coincided  closely  with  the  MKLT  results  so  this  curve 
was  used  for  both.  These  curves  correspond  to  a  1.7  hour  decay  time  when 
the  rate  coefficients  in  question  are  increased  and  to  a  2.7  hour  decay 
time  when  the  coefficients  are  decreased.  The  fact  that  the  effect  on  the 
electron  content  when  the  rate  coefficients  are  decreased  is  less  than  when 
they  are  increased  is  due  to  the  other  reactions,  whose  rate  coefficients 
are  not  varied.  These  become  relatively  more  important  in  the  first  cir¬ 
cumstance  and  limit  the  effect  on  the  overall  time  decay  constant. 
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SECTION  5 

COMMUNICATION  EFFECTS  DUE  TO  CHEMISTRY  VARIATION 


To  determine  the  effects  of  the  chemistry  uncertainties  as  ac¬ 
curately  as  possible  would  involve  running  the  MELT  code  at  least  twice 
with  different  rate  coefficients  and  using  the  results  with  a  communications 
systems  code  for  a  large  number  of  different  design  parameters,  frequencies 

and  transmitter-receiver  configurations.  Not  only  would  this  be  expensive 

2 

but  would  not  resolve  the  effect  of  N(  D)  and  the  lack  of  diffusive  trans¬ 
port  and  neutral  chemistry  in  MELT. 

What  we  have  done  is  to  apply  the  uncertainties  in  the  chemistry 
to  two  values  of  TEC  (tube  electron  content)  which  might  be  encountered 
along  a  transmitter-receiver  path  and  have  tried  to  determine  what  the 
effects  of  these  uncertainties  would  be  on  UHF  and  higher  frequency  com¬ 
munications.  Many  of  these  effects  depend  strongly  on  the  type  of  receiver 
and  its  design  parameters.  We  are  strongly  indebted  to  Robert  Bogusch  for 
enumerating  the  various  communications  effects,  evaluating  a  number  of  these 
effects  using  realistic  frequencies  and  design  characteristics,  and  providing 
the  information  necessary  to  calculate  how  they  would  vary  with  variations 
in  TEC. 


The  values  of  TEC  as  a  function  of  time  used  were  derived  from 
the  MELT  results.  In  Figure  7  is  plotted,  as  the  first  case,  the  TEC  for 
the  central  striation  (upper  solid  line).  We  have  assumed  that  this  approx¬ 
imates  the  TEC  a  ray  path  which  looks  up  the  central  field  line  would  see. 
Actually,  because  of  the  field  line  curvature  this  represents  an  upper  limit 
for  such  a  ray  path.  Ray  paths  probably  could  be  found  (transverse  to  the 
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Fiqure  7.  Tube  electron  content  (TEC)  for  central 

striation  (Case  I)  and  off-axis  striation 
(Case  II). 
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plume  at  low  altitudes),  however,  which  would  exhibit  larger  values  of  TEC. 
The  second  case  was  meant  to  approximate  a  ray  path  looking  up  an  off-axis 
field  line.  For  this  average  TEC  a  curve  like  that  in  the  first  case  was 
used,  but  with  values  uniformly  a  factor  of  three  lower.  This  is  shown  by 
the  lower  solid  line  in  Figure  7.  The  dotted  lines,  together  with  these 
solid  lines,  represent  the  range  of  variation  due  to  chemistry  uncertainties 
and  were  plotted  using  the  decay  times  obtained  in  the  previous  section. 

Before  we  get  into  the  specific  communications  effects  we  can 
examine  the  times  these  various  TEC  curves  take  to  get  into  the  range  of 
values  usually  encountered  in  the  ambient,  natural  atmosphere.  This  is  an 
important  general  consideration  when  dealing  with  any  communication  system.  1 

The  TEC  values  encountered  in  the  natural  atmosphere  are  in  the 

range  of  lO*'5  to  lO1^  cm  “.  We  see  from  Figure  7  that  the  variations  in 

chemistry  cause  significant  uncertainties  in  the  times  the  TEC  takes  to 

decrease  into  the  ambient  regime.  In  the  first  case  (central  tube)  the 

14  -2 

variation  in  time  to  decay  to  10  cm  is  about  4  hours  (from  about  8  1/2 
hours  to  over  12  hours) .  For  the  second  case  (off-axis  tube)  the  variation 
is  almost  3  hours  (from  6  1/2  hours  to  over  9  hours) .  Uncertainties  such 
as  these  could  be  important. 

The  specific  propagation  effects  which  tend  to  degrade  the  per¬ 
formance  of  any  communications  system  can  be  divided  into  two  categories. 

The  first  category  of  effects  include  those  which  depend  only  on  the  mean 
electron  content  along  the  propagation  path.  These  are  absorption,  disper¬ 
sion,  doppler  shift  and  time  delay  rate  of  change.  The  second  category 
includes  scintillation  effects  which  depend  on  the  structure  or  degree  of 
striation  of  the  propagating  medium.  These  are  amplitude  scintillation, 
phase  scintillation,  scintillation  decorrelation  time  and  scintillation 
decorrelation  frequency.  The  estimates  which  follow  regarding  the  magni¬ 
tude  of  these  effects  and  their  variation  due  to  our  chemistry  variations 
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are  crude  estimates  and  depend  on  numerous  simplifying  assumptions,  many 
of  which  will  not  be  explicitly  stated. 


The  absorption,  a,  can  be  expressed  as 
"  r  T 

-  /  n  v  . 

J  e  ei 


a  = 


1.2  x  10 


dZ  dB 


(23) 


where  f  is  the  propagation  frequency  in  I iz ,  Zj  is  the  height  of  the 

transmitter,  and  v  .  is  the  electron-ion  collision  frequency,  v  .  is 
ei  n  1  ei 

given  by 

in  |l.2xl016  sec"1  .  (24) 


v  .  =  1.8  n  T'3/2 
ei  e 


If  we  take  T  =  1000K  and  f  =  300  MHz,  we  obtain 


a  =  1.5  x  10 "22J  n2  dZ  dB 
o 


(25) 


Assuming  the  TEC  is  equivalent  to  a  uniform  n  over  a  layer  thickness 

g  e 

of  10  cm,  we  have 


-30  2 

a  =  1.5  x  io  (TEC) 


dB 


(26) 


If  we  assume  3  dB  attenuation  as  the  criterion  for  significance  we  obtain 
that,  for  f  =  300  MHz,  we  need 

TEC  2  1.4  x  io15  cm"2  .  (27) 

Examination  of  Figure  7  shows  that,  for  the  central  striation  case,  the 
uncertainty  in  chemistry  prolongs  the  time  of  significant  attenuation  from 
about  4  hours  to  5  hours.  For  the  off-axis  striation  case  significant 
attenuation  exists  for  about  2  hours  and  the  chemistry  uncertainty  causes 
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only  about  a  15  minute  effect.  This  last  result  is  strongly  affected  by 
our  assumption  that  the  uncertainty  in  chemistry  starts  at  100  minutes. 

This  was  done  because  of  our  lack  of  specific  knowledge  concerning  the 
various  temperatures  during  the  recompression  phases  of  the  motion  which 
occurs  during  the  first  100  minutes.  There  would,  however,  be  uncertainties 
in  the  chemistry,  primarily  due  to  reaction  19,  during  this  time.  A  time 
of  30  minutes  for  the  prolongation  of  significant  attenuation  would  prob¬ 
ably  be  more  reasonable  for  the  off-axis  striation  case. 

For  L-band  (1.5  GHz j  ,  a  TEC  of  greater  than  7  x  10^  is  needed 
to  produce  significant  attenuation.  This  means  that  attenuation  is  not  a 
problem  at  this  frequency  except  possibly  at  very  early  times. 

Losses  due  to  dispersion  depend  on  the  value  of  the  TEC  and  on 

the  modulation  rate.  At  UHF  (300  MHz),  dispersive  losses  are  negligible 

for  any  TEC  below  10lf)  cm  “  as  long  as  the  modulation  rate  is  less  than 

300  kbps,  a  rate  in  excess  of  any  rate  currently  being  used.  At  L-band 

(1.5  MHz)  and  a  10  Mbps  modulation  rate  the  dispersive  loss  is  about  1  dB 
16  “2 

at  a  TEC  of  10  cm  .  This  is  too  small  a  loss  to  be  significant.  At 

16  2 

X-band  (7.5  GHz)  dispersion  is  negligible  for  any  TEC  below  10  cm  and 
modulation  rates  below  40  Mbps. 

The  degradations  caused  by  doppler  shift  and  time  delay  rate  of 
change  are  also  negligible  for  the  cases  shown  in  Figure  7.  Thus  in  the 
first  category  of  propagation  effects,  only  UHF  absorption  is  important 
and  can  be  significantly  affected  by  our  chemistry  uncertainties  for  the 
cases  we  have  chosen. 

In  the  area  of  scintillation  effects,  more  assumptions  are  needed 
than  in  the  previous  category  of  effects  to  obtain  any  sort  of  quantitative 
results.  One  of  these  assumptions  we  have  made  is  that  the  striations  are 


fully  developed  after  two  hours  and  the  square  of  the  electron  density 
standard  deviation  is  given  by 
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T  n 
4  e 


2.5  x  10_1/  (TEC)2 


(27) 


where  Z  ,  the  total  thickness  of  the  ionized  layer,  was  again  taken  to 
8  ^ 

be  10  cm.  Another  assumption  is  that  the  propagation  is  perpendicular 
to  the  magnetic  field  and  thus  the  striations.  This  was  assumed  to  make 
the  calculations  simpler  although  the  effects  in  this  case  are  less  severe 
than  when  the  propagation  is  parallel  to  the  striations. 


These  simplified  calculations  using  the  TEC  values  in  Figure  7 
yield  the  result  that,  for  UHF,  amplitude  scintillation  is  saturated  at 
the  Rayleigh  fading  limit  for  all  values  of  TEC.  This  is  true  even  at 
ambient  values  and  the  problem  in  this  case  is,  as  it  is  in  the  natural 
ionosphere,  when  and  how  the  striations  decay.  The  problem  probably  does 
not  involve  uncertainties  in  chemistry  in  any  significant  manner.  The 
size  of  the  region  of  saturated  scintillation  will  change,  however,  with 
variations  in  the  TEC  values  as  a  function  of  time. 


For  L-band  (1.5  GHz),  amplitude  scintil lation  is  saturated  for 
6  1/2  to  9  1/2  hours  for  the  smaller  TEC  case  in  Figure  6  depending  on 
the  chemistry  values  used.  For  the  larger  TEC  case  the  spread  in  time 
is  8  1/2  to  over  12  hours.  This  is  the  same  uncertainty  in  the  time  as 
we  obtained  for  the  TEC  values  to  get  into  the  ambient  regime  and  can  be 
considered  significant. 

For  X-band  (7.5  GHz),  the  analogous  times  are  3  to  4  hours  for 
the  lower  TEC  case  and  5  to  7  hours  for  the  larger  TEC  values.  This  also 
can  be  considered  significant.  In  both  the  L-band  and  X-band,  as  in  the 
UHF  case,  the  size  of  the  saturated  region  will  be  sensitive  to  chemistry 


37 


Phase  scintillation  effects  are  directly  proportional  to  the 
value  of  TEC  thus  a  variation  of  a  factor  of  2  or  more  in  TEC  can  translate 
to  the  same  variation  in  the  magnitude  of  phase  scintillation.  This  can 
be  important,  particularly  to  phase  coherent  receivers,  but  quantification 

1  G 

of  system  effects  critically  involves  many  receiver  design  parameters. 

The  scintillation  decorrelation  time  is  the  most  important  signal 
parameter  affecting  receiver  performance  under  saturated  scintillation 
conditions.1  It  varies  approximately  as  the  inverse  of  the  value  of  the 
TEC  and  thus  is  subject  to  a  variation  of  a  factor  of  two  or  more  after 
about  5  hours.  This  can  be  important  to  system  performance  but  as  with 

1  6 

phase  scintillation,  depends  critically  on  the  system  design  parameters. 

The  decorrelation  frequency  provides  a  measure  of  frequency- 
selective  scintillation  effects.  When  this  is  much  less  than  the  signal 
bandwidth,  these  frequency-selective  effects  can  be  quite  important  to 
system  performance.16  This  is  an  area  where  investigation  is  currently 
being  pursued  but  it  can  be  determined  that  the  decorrelation  bandwidth 
varies  approximately  as  the  inverse  of  the  TEC  squared.16  Thus,  the 
variance  of  a  factor  of  2  to  5  in  TEC  after  5  hours  due  to  chemistry  un¬ 
certainty  translates  into  a  factor  of  4  to  25  in  the  decorrelation  band¬ 
width.  This  can  be  quite  important  to  systems  that  use  large  signal  band- 
widths  (such  as  GPS  and  DSCS  spread  spectrum  links).16 

At  times  greater  than  5  hours,  when  the  variation  due  to  uncer¬ 
tain  chemistry  is  greatest,  X-band  and  higher  frequency  systems  tend  not 
to  be  scintillation  saturated.  Thus  the  sensitivity  to  chemistry  of 
scintillation  effects  is  more  important  for  L-band  and  UHF  systems. 

To  summarize  the  main  areas  where  reasonable  chemistry  variations 
will  affect  communications,  we  have  first  that  the  resultant  variation  in 
times  to  reach  ambient  ionospheric  conditions  is  significant.  The  times 
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for  which  signal  attenuation  due  to  ion-electron  collisions  is  important 
(>  3  dB)  will  also  vary  significantly  for  UHF  systems.  The  same  is  true 
for  the  variation  in  times  over  which  we  have  amplitude  scintillation 
saturation  at  L-band  and  X-band  frequencies.  The  size  of  the  region  of 
scintillation  saturation  will  probably  also  vary  significantly  with  chem¬ 
istry  but  quantification  is  difficult  without  much  more  analysis.  Finally, 
during  the  time  of  scintillation  saturation,  important  degradation  effects 
associated  with  phase  scintillation  and  scintillation  decorrelation  time 
and  frequency  will  probably  vary  significantly  with  chemistry  variations, 
particularly  for  UHF  and  L-band  systems  which  are  saturated  to  late  times. 
These  effects  are  critically  dependent,  however,  on  specific  systems  and 
receiver  design  parameters. 
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SECTION  6 

RECOMMENDED  RATE  COEFFICIENTS 


The  most  important  reaction,  reaction  19,  is  the  most  difficult 
one  for  which  to  determine  a  coefficient.  The  problem,  as  indicated  before, 
is  that  there  is  no  experimental  data  available  when  both  the  translational 
and  nitrogen  vibrational  temperatures  are  above  about  luOOK.  This  temperature 
region,  however,  is  where  most  of  the  action  takes  place  following  high 
altitude  bursts.  We  have  constructed  a  rate  coefficient  based  on  several 
assumptions.  The  first  of  these  is  that  the  rate  coefficients  and  transla¬ 
tional  temperature  dependence  of  the  nitrogen  V  =  0  and  V  =  1  vibrational 
states  are  identical.  This  probably  isn't  true,  but  is  consistent  with  the 
fact  that,  at  300K  translational  temperature,  the  rate  constant  does  not 
increase  with  vibrational  temperature  until  the  V  =  2  state  is  significantly 
populated.  The  second  assumption  is  that  the  rate  coefficients  for  the 
higher  vibrational  states  are  independent  of  translational  temperature. 

These  are  taken  to  be  5  x  10  1.2  x  10  2.5  x  10  ^  and  4.0  x  lcf 10 

3  12 
cm  /sec,  respectively,  for  the  v  =  2,  v  =  3,  v  =  4,  and  v  >  4  states. 

The  assumption  of  no  variation  with  kinetic  temperature  is  probably  not 

true  either  but  seems  to  be  the  most  reasonable  assumption  in  the  absence 

of  any  hard  experimental  data.  The  third  assumption  is  that  the  vibrational 

levels  in  nitrogen  are  equally  spaced.  In  view  of  the  uncertainty  caused 

by  the  first  two  assumptions,  the  error  involved  in  this  is  trivial.  Using 

these  assumptions  and  some  curve  fitting  we  obtain 

/  -6600/Ty\  -6600/Tv 

kin  =  k  ll-e  )+  k  e  cm'Vsec  (23) 

19  tr  V  /  v  v 
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where 


k  =  1.2  x  10 
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and 


T  \0.8 


k  =  9.6  x  10 
v 


n(M 

\3000 / 


(25) 


and  T  and  T  are  the  vibrational  and  translational  temperatures, 
v  tr 

For  reaction  16,  the  N+  +  0  charge  exchange  reaction,  we  have 
already  indicated  our  guess,  i.e.. 


k,.  =  3  x  io 

lo 


■13 


(m) 


3, 

cm  /sec 


(26) 


where  T  is  the  effective  translational  temperature. 


★ 

More  intormation  is  now  available  on  this  reaction  rate  coefficient.  The 
analysis  of  the  Atmospheric  Explorer  satellite  data  has  yielded  a  value 
of  2.2  x  io* 13  cm'Vsec  at  840K  with  a  claimed  probable  accuracy  of  ±20  per¬ 
cent.  This  is  significantly  larger  than  the  value  obtained  by  extrapolating 
the  high  energy  data.  It  now  appears  this  rate  coefficient  may  behave  as 
many  of  the  other  ion-molecule  coefficients,  i.e.,  first  decreasing  with 
increasing  temperature  and  then  increasing  after  going  through  a  minimum. 

We  would,  therefore,  now  recommend  (see  Figure  1) 


_  -12 

/300 1 

1 1/2  3  . 

16  =  3.69  x  io 

l  T  j 

'  cm  /sec 

T 

<  1600K 

16  =  3.0  x  io’13  ( 

'  T  \ 
300  1 

3/ 

cm  /sec 

T 

>  1600K 

This  is  an  increase  of  a  factor  of  2  at  1000K  over  the  value  we  have  used. 
Since  we  used  a  rather  conservative  uncertainty  in  our  calculations  and 
since,  as  we  have  pointed  out,  an  increase  in  this  rate  would  not  have  much 
of  an  effect  at  times  after  stabilization,  our  results  and  conclusions  are 
not  seriously  affected  by  this  new  data. 
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For  reaction  18,  the  0  +  0^  charge  exchange  reaction,  we  have 

examined  the  data8’9  ’ 1 0  ,n’ 1 4  and  obtained 


k18  '  2'‘  X  10 


k18  -  7.6  x  10 


f-T)' 

iiof 


cm  /sec  T  <  1800K 


cm  /sec  T  >  1800K 


where,  as  before,  T  is  the  effective  translational  temperature. 


For  the  sum  of  reactions  14  and  15  a  value  of  6  x  io~  from 

300K  to  6000K  is  indicated  by  the  data.9’10’11  Above  6000K  the  rate  coef- 

1/2 

ficient  increases  with  a  T  '  dependence. 10  At  these  high  temperatures, 
however,  there  is  little  likelihood  of  being  present  in  any  appreciable 

amount  so  the  high  temperature  dependence  is  unimportant. 

We  have  discussed  the  possibility  that  NO  may  be  present  in 
the  plume  in  significant  quantities.  Unfortunately,  the  MELT  code  does 
not  carry  any  of  the  neutral  reactions  and  processes  (e.g.,  diffusion) 
involved  so  it  is  difficult  to  even  estimate  the  NO  production.  If, 
however,  appreciable  amounts  of  NO  can  be  found,  certain  other  ion- 
neutral  reactions  should  be  considered  for  their  possible  importance. 

Two  of  these  reactions  are 


0+  +  NO  ->  N0+  +  0 


N+  +  NO  -*•  N0+  +  N 


The  reaction  28  rate  coefficient  has  a  value  at  300K  of  about  8  x  10'13 
and  behaves  at  higher  temperatures  in  somewhat  a  similar  manner  as  the 
reaction  19  rate  coefficient,  although  it  does  not  dip  down  as  strongly 
between  300K  and  1200K  nor  rise  as  sharply  above  1200K.  Since  it  is  highly 


unlikely  that  the  NO  concentration  will  ever  be  comparable  with  the  N-, 

concentration,  this  reaction  is  probably  unimportant.  Reaction  29,  however, 

.-9  3 

has  a  rather  large  rate  coefficient.  At  300K  the  value  is  about  10  cm  /sec 
and  it  is  unlikely  that  it  changes  significantly  above  300K.  Thus  if  the 

_3 

NO  concentration  is  of  the  order  of  10  of  the  atomic  oxygen  concentra¬ 
tion  it  can  compete  with  reaction  16  at  ambient  temperatures  and  directly 
convert  N+  to  a  molecular  ion. 

Another  reaction  which  must  be  considered  is 

0+  +  N(2D)  -*  0  +  N+  .  (30) 

Analysis  of  the  Atmospheric  Explorer  satellite  data17  has  yielded  a  value 

of  1.3  *  lO1^  for  the  rate  coefficient  of  reaction  30  and  it  is  claimed 

that  this  reaction  is  important  in  the  natural  ionosphere. 17  What  the 

effect  would  be  in  the  nuclear  disturbed  case  is  difficult  to  ascertain. 

2 

If  the  N(  D)  concentration  is  within  two  orders  of  magnitude  of  the  N2 
concentration,  the  0+  will  be  converted  more  rapidly  to  N+  by  reaction  30 
than  to  N0+  by  reaction  19. 
and  to  make  reactions  14,  15, 

0)  relatively  more  important 


This  would  tend  to  slow  the  total  ion  decay 
and  16  (the  reactions  of  N+  with  0^  and 
than  we  have  assumed. 
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SECTION  7 

RECOMMENDATIONS  FOR  FURTHER  STUDY 


The  obvious  recommendation  would  be  to  reduce  the  rate  coefficient 
uncertainties  which  have  been  determined  to  be  important.  The  problem  is 
that  one  of  the  reasons  these  uncertainties  still  persist  long  after  their 
importance  was  recognized  is  the  extreme  difficulty  in  doing  experiments 
under  the  conditions  of  interest.  Reaction  16,  the  N+  +  0  charge  exchange 
reaction,  which  is  important  at  all  altitudes  within  the  plume  below  about 
500  km,  is  a  case  in  point.  The  interesting  region  is  from  0.5  eV  energy 
down  to  thermal  energies.  The  measurements,  however,  do  not  go  below  0.5  eV 
and  these  were  not  easy  to  obtain  and  have  been  accomplished  only  recently. 
Fortunately,  the  decrease  in  the  rate  coefficient  of  reaction  19  at  ambient 
conditions  below  what  was  used  in  the  MELT  code  makes  it  more  improbable 
that  reaction  16  will  control  the  total  ion  decay  at  any  time  after  plume 
stabilization.  A  decrease  of  two  or  three  in  the  rate  coefficient  of  the 
N+  +  0  reaction  at  1000K  from  the  value  obtained  by  extrapolating  the 
experimental  values  would  be  needed  to  make  this  reaction  the  controlling 
factor.  Since  the  extrapolation  involves  a  range  in  temperature  of  a  factor 
of  4  rather  than  the  factor  of  40  we  had  prior  to  the  recent  measurements, 
this  large  a  decrease  in  the  rate  coefficient,  though  possible,  is  not  very 
likely.  At  times  prior  to  stabilization,  when  the  ion  temperature  may  be  as 
high  as  1  eV,  the  0+  +  reaction  goes  much  faster  than  at  ambient  temper¬ 
ature  and  the  N  +0  reaction  will  probably  be  the  rate  controlling  re¬ 
action,  particularly  since  the  0^  is  depleted  and  diffusion  certainly 
has  not  had  time  to  replace  it.  At  these  higher  energies  we  are,  however, 
in  or  very  near  the  range  of  the  experimental  data.  The  uncertainty, 
therefore,  is  less  than  in  the  region  of  the  data  extrapolation.  It  would 
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be  reassuring  to  have  experimental  data  in  the  region  of  0.1  and  0.5,  par¬ 
ticularly  at  the  low  end  of  this  energy  range,*  but  this  is  not  of  the 
highest  priority.  The  uncertainty  in  reaction  18,  the  0+  +  0.,  charge  ex¬ 
change  reaction,  which  is  important  below  about  350  km  altitude,  has  been 
shown  to  have  some  measure  of  importance.  This  importance  is  not  as  great 
as  the  importance  in  the  uncertainty  in  reaction  19  except  possibly  at  very 
late  times.  In  addition,  this  uncertainty  stems  largely  from  the  scatter 
in  the  data  from  a  number  of  experiments.  Assuming  all  the  experiments 
are  equally  reliable,  which  they  are  not,  a  reduction  in  this  uncertainty 
would  involve  doing  a  large  number  of  further  experiments.  Finally,  this 
uncertainty  is  probably  overshadowed  by  uncertainties  introduced  by  the 
lack  of  0-,  diffusion  and  neutral  chemistry  in  the  plume.  For  these 
reasons  reducing  the  uncertainties  in  the  rate  coefficient  of  reaction  18 
is  also  not  of  the  highest  priority. 

The  uncertainty  in  the  rate  coefficient  of  reaction  19,  which  is 
important  at  all  altitudes  up  to  500  km,  is  the  most  important  uncertainty 
we  have  discussed.  This  is  the  rate  controlling  reaction  for  total  ion 
decay  and  occurs  at  temperatures  not  encompassed  by  the  experimental  data. 

It  is  of  high  priority  to  measure  this  rate  coefficient  when  both  kinetic 
and  N'2  vibrational  temperatures  are  above  1000K.  In  particular,  measure¬ 
ments  which  would  reveal  the  dependence  on  kinetic  temperature  of  the  rate 
coefficient  for  the  reaction  when  N2  is  in  the  V  =  1,  2,  and  3  vibra¬ 
tional  states  would  be  extremely  useful. 

2 

The  uncertainty  of  N(  D)  depends  strongly  on  the  uncertainties 
in  the  excited  state  production  in  the  initial  energy  deposition  phase, 
the  production  by  electron  impact  during  the  recompression  phases  of  plume 
development,  and  the  quenching  by  atomic  oxygen.  The  quenching  is  particularly 

*  -12  3 

The  AE  satellite  data  analysis  has  yielded  a  value  of  2.2  x  10  cm  /sec 

at  840K. 17 


important  and  of  high  priority.  If  the  quenching  rate  coefficient  is  of 
-11  2 

the  order  of  10  ,  the  effect  of  N(  D)  will  be  restricted  to  only  the 

first  minute  or  two  after  its  production,  either  initially  or  during  plume 

recompression,  and  thus  not  very  important  to  communications.  If  the  rate 

- 12  2 

coefficient  is  of  the  order  of  10  or  less,  the  effect  of  N(  D)  might 
be  important  to  communications  for  times  of  the  order  of  an  hour  after 
production.  A  measurement  of  this  rate  coefficient  is  therefore  of  high 
priority. 


2 

In  this  respect,  subsidiary  questions  arise.  When  N(  D)  is 
quenched  by  atomic  oxygen,  is  an  0(^D)  state  produced?  Also  when  NO 
dissociatively  recombines,  is  an  Of^D)  produced  in  the  channel  where  the 

4 

nitrogen  is  produced  in  the  N(  S)  state?  If  the  answers  to  both  of  these 
are  yes,  then  we  may  have  to  worry  about  the  reaction 

N(4S)  +  0(1D)  ->  N0+  +  e  .  (31) 

Quenching  of  0(^D)  by  N_,  which  has  a  rate  coefficient  of  about 

1  1  *7  ^ 

3  x  10-1  cni  /sec,4  as  well  as  a  small  rate  coefficient  for  reaction  31, 

probably  makes  this  unimportant.  However,  a  large  production  of  0(^D) , 

and  its  subsequent  quenching  by  N^,  which  results  in  N2  vibrational 

excitation,  could  substantially  affect  the  N2  vibrational  distribution. 

This  would  impact  strongly  on  the  value  of  the  rate  coefficient  of  reaction  19. 

Further  investigation  is  needed,  but  it  would  be  useful  if  the  amount  of 
1  2 

production  of  0(  D)  for  the  N (  D)  quenching  reaction  were  measured  as 
well  as  the  value  of  the  rate  coefficient. 


A  further  possible  effect  of  N(2D)  stems  from  the  reaction 
0+  +  N(2D)  ->•  N+  +  0  ( 

which,  according  to  the  Explorer  satellite  data  analysis,  has  a  fairly 
large  rate  coefficient.  The  effect  of  this  reaction  on  the  total  ion 
decay  should  be  investigated. 
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The  last  uncertainty  we  have  discussed  is  that  involved  in  the 
depletion  of  0^  in  and  near  the  plume,  its  replacement  by  diffusion, 
and  the  mitigation  of  its  replacement  due  to  neutral  reactions  involving 
atomic  nitrogen.  This  impacts  not  only  on  all  ion  reactions  involving  0^ 
but  on  the  production  of  NO,  whose  presence  and  persistence  could  sub¬ 
stantially  alter  the  ambient  ionosphere  for  very  long  periods  of  time  after 
a  number  of  high  altitude  bursts.  This  could  involve  important  communica¬ 
tions  effects  and  should  be  investigated. 
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